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ABSTRACT 


The  development  of  a  speed-of-light  hard-kill  weapon  system  for  military 
applications  represents  a  significant  advancement  in  technology  over  present 
conventional  kinetic  weapon  systems.  Over  the  past  two  decades,  the  US  Navy  has 
successfully  developed  a  megawatt-class  chemical  laser;  however,  under  some  maritime 
environments,  the  high  power  beam  propagation  was  unable  to  delivery  sufficient  energy 
to  kill  a  modern  anti-ship  missile  (ASM)  due  to  significant  atmospheric  absorption  and 
the  resulting  thermal  blooming  process.  A  critical  problem  to  resolve  for  the  shipboard 
high-energy  laser  weapon  systems  is  to  develop  a  shipboard-compatible  megawatt-class 
laser  weapon  at  a  wavelength  where  the  atmospheric  absorption  is  smallest.  The 
megawatt-class  Free  Electron  Laser  (FEL)  has  significant  advantages  over  conventional 
weapon  systems  and  other  chemical  high-energy  laser  systems.  Infinite  magazine,  rapid 
response,  and  wavelength  tunability  make  the  FEL  a  suitable  and  desirable  shipboard 
weapon  system. 

This  thesis  divides  into  four  chapters.  Chapters  I  and  II  introduce  the  FEL  and 
background  theory  of  the  FEL.  Chapter  III  explores  the  analysis  of  the  LANL 
Regenerative  MW  FEL  Amplifier  design  and  optimizes  its  efficiency.  Lastly,  chapter  IV 
summarizes  the  feasibility  of  achieving  the  desired  efficiency. 
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I.  INTRODUCTION 

A.         DEFENSE  AT  THE  SPEED  OF  LIGHT 

In  today's  technology,  weapon  systems  rely  predominantly  on  kinetic  kill,  striking 
a  missile  with  another  missile  or  with  several  bullets.  These  weapons  are  relatively  slow 
and  require  a  longer  detect-to-engage  sequence  when  factoring  in  tracking  time,  launch 
time  and  time  of  flight  to  intercept.  In  searching  for  faster  and  faster  weapons  and  ever- 
decreasing  response  times,  Directed  Energy  Weapons  (DEW)  have  been  on  the  U.S. 
military's  "love  to  have"  list  for  several  decades.  DEWs  are  attractive  for  several  reasons, 
and  the  most  predominant  reason  is  that  they  offer  speed-of-light  defense. 

In  1971,  the  U.S.  Navy  established  the  Navy  High  Energy  Laser  (HEL)  Program 
with  the  intentions  of  developing  speed-of-light  hard-kill  weapon  system.  Initial 
development  focused  on  building  a  CO2  gas  dynamic  laser  technology,  but  a  few  years 
later  its  emphasis  shifted  to  continuous-wave  deuterium  fluoride  (DF)  chemical  lasers 
because  they  propagate  far  better  in  the  maritime  environments.  In  the  early  1980s,  the 
megawatt  class  Mid-Infrared  Advanced  Chemical  Laser  (MIRACL)  was  built  by  TRW 
for  the  Navy  as  the  highest  average  power  continuous  wave  (CW)  laser  in  the  west.  The 
SeaLite  Beam  Director  (SLBD)  was  built  by  Hughes  Aircraft  to  accommodate  and  direct 
the  MIRACL  beam.  The  MIRACL  and  SLBD  have  successfully  tested  at  White  Sands 
Missile  Range,  New  Mexico,  in  engaging  and  destroying  several  types  of  missiles  in 
flight  [1]. 

Although  the  MIRACLE/SLBD  system  was  successfully  tested  against  crossing 
targets  in  a  fleet  defense  scenario,  its  lethality  was  not  sufficient  in  some  ship  self-defense 
engagement  scenarios  due  to  insufficient  crosswind  [1].  Insufficient  crosswind  will  cause 
the  laser  beam  heating  the  air  in  its  own  path  resulting  in  a  defocusing  optical  lens  known 
as  thermal  blooming  [2J.  The  thermal  blooming  effect  can  severely  reduce  the  High 
Energy  Laser  Weapon  Systems'  (HELWS)  lethal  range.  As  a  result,  it's  necessary  to 
replace  the   MIRACL  with  a  megawatt-class  laser  device  at  a  wavelength  where 


atmospheric  absorption  is  smallest.  An  ideal  solution  is  the  megawatt-class  Free  Electron 
Laser  (FEL)  that  is  tunable  to  the  optimum  propagation  wavelength. 

B.         SHORT  HISTORY  OF  FREE  ELECTRON  LASER 

The  first  generation  of  coherent  radiation  from  free  electron  beams  was 
microwave  tubes  that  relied  on  slow-wave  structures,  and  their  development  received  a 
significant  impact  from  the  radar  development  during  World  War  II.  In  1951,  Motz 
showed  that  an  electron  beam  propagating  through  an  undulator  magnet  can  be  used  to 
amplify  radiation.  After  1 960,  research  on  short  wavelength  lasers  has  been  dominant,  but 
they  require  an  atomic  or  molecular  medium  to  operate,  resulting  in  a  limitation  on 
wavelength  tunability.  In  1971,  J.  M.  J.  Madey  of  Stanford  University  brought  attention  to 
what  he  proposed  as  the  'free  electron  laser'.  Relying  on  research  of  synchrotron  radiation 
sources,  Madey  conceived  a  device  that  would  lase  in  the  visible  using  a  beam  of 
relativistic  electrons  as  the  source.  In  1976,  the  first  successful  free  electron  laser  (FEL) 
experiment  was  demonstrated  by  Madey  and  his  co-workers  at  Stanford  University,  in 
that  they  measured  7%  gain  from  an  FEL  configured  as  an  amplifier  at  10 jam 
wavelength.  One  year  later,  another  successful  FEL  operation  of  the  same  FEL  configured 
as  an  oscillator.  These  successful  experiments  created  a  large  interest  in  FEL  research. 
Since  then,  scientist  and  experimental  groups  around  the  world  built  FELs  that  operate  at 
frequencies  ranging  from  microwave  to  the  UV  [11].  Three  other  FELs  were  built  soon 
after.  The  first  was  built  at  Laboratoire  pour  l'Utilisation  du  Rayonment 
Electromagnetique  in  Orsay,  France,  which  operated  in  the  visible  range.  The  second 
involved  a  team  from  TRW  at  Stanford  using  the  Superconducting  Accelerator  (SCA), 
lasing  in  the  near  infrared  (IR).  The  third  was  built  at  Los  Alamos  that  lased  in  the  mid- 
infrared  [3].  Now,  there  are  many  FELs  used  for  research  around  the  world. 


c. 


BASIC  FREE  ELECTRON  LASER  THEORY 


The  FEL  is  a  device  that  transforms  the  kinetic  energy  of  a  relativistic  electron 
beam  into  electromagnetic  (EM)  radiation  [4].  This  is  accomplished  by  sending  a 
relativistic  electron  beam  from  a  linear  accelerator,  microtron  or  a  storage  ring  etc., 
through  an  alternating  magnetic  field  produced  by  a  device  known  as  an  undulator  as 
shown  in  Figure  1.1.  The  interaction  of  the  negatively  charged  electrons  with  the  spatially 
periodic  magnetic  field  induces  a  transverse  acceleration  on  the  electrons  that  causes 
them  to  exchange  energy  with  a  copropagating  radiation  field.  The  kinetic  energy 
extracted  from  the  electron  beam  is  converted  to  electromagnetic  radiation.  The  energy  of 
the  electron  beam  and  the  period  of  the  undulator  magnet  determine  the  wavelength  of 
radiation.  This  system  provides  a  powerful  source  of  tunable,  coherent,  electromagnetic 
radiation  that  can  operate  in  a  wavelength  range  from  about  1  cm  to  the  X-ray  region. 
This  system  can  also  provide  continuous  wave  (CW)  operation,  or  very  short  pulses  in  the 
nanosecond,  picosecond,  or  even  subpicosecond  region. 


Electron  Beam 


Undulator 


Resonator 


Mirror 


Figure  1.1.  Schematic  diagram  of  a  free  electron  laser. 


II.  PRINCIPLES  OF  OPERATION  OF  FEL 

A.         ELECTRON  TRAJECTORY  IN  THE  UNDULATOR 

For  better  understanding  of  the  physical  mechanism  underlying  the  operation  of  an 
FEL,  we  first  consider  the  interaction  of  a  free  electron  in  the  undulator  and  an 
electromagnetic  wave.  The  Lorentz  force  equation  and  the  equation  for  the  energy 
exchange  between  the  electron  and  the  electromagnetic  wave  are 


dp 
dt 


'-     v      -^ 

E+-xB 

c 


(2.1) 


d(ymc2) 

where  e,  m  and  v  are  the  charge,  mass  and  velocity  of  the  electron,  respectively;  y  is  the 
relativistic  factor  of  the  electron,  y-  l/^/(l-/?2  J ;  /?  =  v/c;  p  is  the  relativistic 
momentum  of  the  electron,  p=ymv  ;  E  is  the  electric  field  of  the  EM  wave;  B  is  the 
magnetic  field  ;  ymc2  is  the  electron  energy;  and  c  is  the  speed  of  light  in  vacuum  [5]. 

Considering  an  electron  with  velocity  v  =  ft  c  and  a  transverse  electromagnetic 
wave  propagating  colinearly,  in  the  z-direction  in  a  vacuum,  then  equations  (2.1)  and 
(2.2)  become 

^P-  =  -—(E  +  pxB),  (2.3) 

dt  mc 

dy         e     -    _ 

-j-=  —  p-E.  (2.4) 

dt        mc 

Assume  the  undulator  has  a  helical  magnetic  field  of  the  form  B=  B[cos(k0z),sin(k0z),Q], 
where  k0  -  2n  I  X0  is  the  undulator  wavenumber,  and  A0  is  the  undulator  wavelength. 
The  corresponding  circular-polarized,  plane-wave  optical  field  is  then 

Es  =  E  (cos if/,-  sin  y,0),  and  Bs=E(sin y/,cos y/,0) , 
where  y/  =  kz  -  cot  +  <f>  .  E  is  the  electric  and  magnetic  field  amplitude  in  cgs  units,  t  is 
time,  (j)  is  the  optical  phase,  k-colc  is  the  optical  wavenumber,  and  (O  is  the  optical 
angular  frequency.  Substituting  the  fields  into  Lorentz  force  equations  yields 
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d(yB.)         e 

-^7T-=— [£(l-/U(cos^,-siny/,0)  +  p,B(-sm(kozUos(koz),0)l    (2.5) 
at  mc 

— ^  =  —  [E(fixcosy-fix  siny/)  +  £(/?>n(£0z)-/?rcos(*0z))],    (2.6) 

and       -— =- —  E(Bxcosu/-Bvsmw),  (2.7) 

a/        mc 

where  /?  =  /?±  +  /?.,   /?±  =  /?rz  +  ft  j ,  and   /?,  =  /?_&.  For  relativistic  electrons,  /?,  «1, 

hence  £(1  -  >9,)  «  ft,  B,  so  in  the  right-hand  side  of  (2.5)  the  first  term  is  ignored.  The 
second  term  is 

d(YpL  )=_  eBf^_  (_ sin(koZ)cos(koZ)Q)  .  (2.8) 

dt  mc 

Integrating  the  transverse  equation  (2.8)  and  assuming  a  perfect  injection  of  electron 

beam  into  the  undulator  cavity,  which  makes  the  constant  of  integration  zero,  yields 

yP±=-\-^T)(cos(koz),sin(koz),0) .  (2.9) 

\mc   J 

Divide  both  side  of  (2.9)  by  /  to  get 

0jz)=  --(cos(koz),sm(koz),0),  (2.10) 

r 

where      the   undulator  parameter  is    K=eB  I  k0mc2  =eBA0  1 2/zmc2  .   Equation  (2.10) 
describes  the  transverse  motion  of  electrons  in  the  helical  undulator  [6]. 


B.    FEL  RESONANCE  «  ELECTRON  AND  PHOTON  RACE 

Accelerated  electrons  traveling  at  a  relativistic  speed  along  the  longitudinal  axis  of 
the  undulator  of  the  FEL  radiate  optical  energy  [6].  In  order  for  the  FEL  to  have  gain  (i.e. 
the  fractional  change  in  optical  power  in  a  single  pass  through  the  undulator),  a  net 
transfer  of  energy  from  the  electron  beam  to  the  optical  beam  must  occur.  In  order  for  an 
optimum  exchange  of  energy  from  an  electron  to  the  optical  field  to  occur,  the  electron 
must  oscillate  in  phase  with  the  optical  field.  This  optimum  energy  exchange  occurs  when 
the  electron  falls  behind  the  optical  fields  by  a  distance  of  one  optical  wavelength  within 
the  span  of  one  undulator  period  [6]. 
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Figure  2.1.  Electron-Photon  race. 


Consider  an  optical  wave  that  travels  at  speed  c  and  passes  over  an  electron  traveling  at 
speed  Pzc  as  shown  in  Figure  2.1.  At  resonance  condition,  one  wavelength  of  the  optical 
wave  passes  over  an  electron  while  it  travels  through  one  undulator  wavelength  X0 ,  so 

that  light  wavelength  X=  (c  -  fizc)t .  Insert  t=  XJ (5zc  yielding  X=—- —  X0 .  Substitute 


A 


into  y~2  =  1-  p2  =  1-  p]  -  Pi  yields 


«v 


(\  +  K2)y-2=\-pz2  . 


(2.11) 


\  +  K2 
Since  /?,  »1,  I-/?,  ~  ~ — — ,  then 

\  +  K2 

2y 

The  equation  (2.12)  is  the  resonance  condition  that  demonstrates  the  tunability  of  the 
FEL.  By  adjusting  the  undulator  parameter  K,  or  the  energy  of  the  electron  beam 
(y  -\)mc2  or  the  undulator  wavelength  A0,  the  wavelength  of  the  optical  field  can  be 
tuned  to  a  desired  wavelength  within  a  prescribed  bandwidth  [7].  This  is  a  significant 
advantage  over  conventional  lasers  that  are  typically  constrained  to  a  fixed  wavelength. 

C.         PENDULUM  EQUATION 

This  section  derives  the  equation  of  motion  of  the  electrons  resulting  from  the 
interaction  with  the  optical  field.  Substitute  (2.10)  into  (2.7)  to  get 

/=— cos(£  +  0),  (2.13) 

ymc 

where 

£  +  </>  =  y/  +  k0z  =  {k  +  k0)z-G)t  +  <j) .  (2.14) 

For  -;r/2<<£'  +  ^<;r/2,the  cosine  function  in  (2.13)  is  positive  and  the  electrons  will 
gain  energy  from  the  optical  field,  whereas  for  n  1 2  <  £  +  <fi  <  3k  1 2 ,  the  cosine  is 
negative  and  the  electrons  will  lose  energy  to  the  optical  field.  Differentiating  (2.1 1)  with 
respect  to  time  t  gives  2y~*y{[  +  K2  )=  2/3 :  ft, .  Relating  y  to  0Z  gives 

t=zlAK,  (2.i5) 

y     (\  +  K2) 
Differentiating      (2.14)      once      yields       C,  =  (k  +  k0)z-co  =  (k  +  k0)c/3,  -co,      and 
differentiating  twice  yields  £  =  (k  +  k0)cfi, .  Insert  the  resulting  expression  for/?,   into 
(2.15)  to  get 

f  r2P£ 


y      (\  +  K2)(k  +  k0)c 


(2.16) 


Since  A  «  A0,  then  k0  «  k,  so  (k  +  k0)c*  kc  =  w.  Using  this  approximation  and  the 
resonance  condition  (2.12),  equation  (2.16)  becomes 

^~ — ^4 — =—  ■  (2-17) 

y      (\  +  K  )co     2co0 

£       Y       eKE 
So,  - —  =  —  =  —z — cos(  C,  +  <f>),  then  the  electron  equation  of  motion  is 
2coo      y      Y'mc 

2coneKE 

£  =  — \ cos(^  +  (f>).  (2.18) 

Y~mc 

Equation  (2.18)  describes  the  electron  phase  dynamics  in  the  form  of  the  pendulum 

equation  [6]. 

Consider  time  evolution  through  the  undulator  L/  f30c  for  an  average  electron, 

where   L   is  the   length  of  the   undulator.    Introducing  relevant  dimensionless  time 
r  =  P0ct  I L  «  ct  I L  ,  we  see  that  r  will  vary  between  0  and  1  along  the  undulator  length. 

Substituting  dt  =  (LI c)dz  into  the  pendulum  equation  (2.18)  gives 

^  =  v  =  |a|cos(^  +  ^),  (2.19) 

where  the  notation  (..)  =  d(..)l dr  ,  \a\  -  AnNeKLE I Y2mc2\s  the  dimensionless  optical 
field  strength,  and  N  =  LI  A0  is  the  total  number  of  undulator  periods.  Equation  (2.19)  is 

the  electron  pendulum  equation  that  governs  the  phase-space  motion  of  the  electrons 
within  the  undulator  and  under  the  influence  of  the  optical  wave  [6]. 
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D.         OPTICAL  WAVE  EQUATION 

An  electromagnetic  wave  will  propagate  in  the  absence  of  a  current  density 
according  to 


V2- 


1    d 


2  N\ 


A(r,t)  =  0,  (2.20) 


J 


c2  dt2 

where  the  Laplacian  is  V2  =  d]  +  d2 .  +  d\ ,  r  is  the  vector  position,  t  is  time,  and  A  is  the 

optical  vector  potential  [5].  The  optical  electric  and  magnetic  fields  can  be  expressed  in 
terms  of  the  vector  potential  by  the  following  relationships 

E  =  -~,  (2.21) 

c  dt 

and       B  =  VxA.  (2.22) 

Assuming  that  the  optical  beam  is  a  circularly  polarized  plane  wave,  then  the  vector 

potential  has  the  form 

A  =£(^)(sinv|/,cosM/,0)=^g  j-3,(r,0&?'(fe"ffl0  1 ,  (2.23) 

where  (£(r,t)  =  E(F,t)e'4'(r'') ,  y/  =  kz  -  cot  +  <j>  is  the  plane-wave  phase,  e  =  (-z',1,0)  is  the 

polarization  vector,  k  =  2ttc  I  co  is  the  optical  wave  number,  co  is  the  optical  frequency, 
and  </>(r,t)  is  the  optical  phase  [7].  For  convenience,  we  write 

A  =  -&(r,t)ee'a  ,  (2.24) 

k 

where  a  =  kz  -cot .  Differentiating  (2.24)  with  respect  to  t  gives 

^  =  -(<r-u»<rVa,  (2.25) 

dt      kK  r 

where  (..)=d(..)  /  dt .  Differentiating  (2.25)  with  respect  to  /  gives 


d2A      e 


=  -(^'-2ico^-co2^ya,  (2.26) 


dt2       k 
where   (..)=d2(..) /  dt2 .  Differentiate  (2.24)  with  respect  to  z  gives 
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— =i(*"-i**>ia, 

dz      k 
where  (..)'=d(..)/  dz  .  Differentiate  (2.27)  with  respect  to  z  gives 

^  =  -(<r"+2ik<r'-£2<ry<\ 

dz2     ky  r 


(2.27) 


(2.28) 


where  (..)"  =d2  (..)/  dz2 .  Since  the  laser  beam  consists  of  a  narrow  range  of  frequencies, 
the  variation  in  the  electric  field  amplitude  and  phase  with  respect  to  space  over  an 
optical  wavelength  X  is  small,  so  that  <f '«  £<f  and  $'  «k<f> .  The  variation  in  the  electric 
field  amplitude  and  phase  with  respect  to  time  over  an  optical  period  is  also  small,  so  that 
<9  «co  <9      and  (f>«a>0    [7].  Therefore,   J"' is  negligible  compare  to  k$' ,  and  J3    is 

negligible  compare  to  co<f  ,  then  (2.26)  can  be  simplified  to 


d2A     ee'c 


dr 


{-li&it-eo2^), 


(2.29) 


and  (2.28)  can  be  simplified  as 


d2A     ee' 


(2ik<T-£2<r). 


dz2        k 
Substitute  (2.29)  and  (2.30)  into  the  left  hand  side  of  (2.20)  yields 


(2.30) 


i  a 


2  A 


V- 


c2  dt 


where  V2,  = — --\ .  Hence, 

1    ax2    dy2 


A{r,t)  = 


ee 


V2i&'  +  2ik\g"  +  -& 
V         c     J 


(2.31) 


{dz   cdt. 


ar  =  o. 


(2.32) 


The  first  partial  second  derivatives  with  respect  to  x  and  v  describe  how  the  beam  changes 
due  to  diffraction  as  it  propagates.  The  last  partial  derivatives  with  respect  to  z  and  t 
describe  the  propagation  of  a  plane  wave  along  the  z  axis. 

Introduce  dimensionless  time  r  =  ct  I L  so  that  r  varies  from  0  to  1  along  the 
undulator  length.  It's  convenient  to  follow  the  wavefront  with  a  new  coordinate 
z  =  z-ct .  With  this  coordinate  change, 
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'dz}  d 


dz     {dz 
since  dz  /  dz  =  1  and  dr  /  dz  -  0  .  Also, 


dz 


+ 


Kdz  j 


dr      dz 


(2.33) 


c  dt      c 


fdz) 
ydtj 


di 


+ 


Kdt  j 


dr 


_d_  j__a_ 

dz      L  dr 


(2.34) 


since   dz  /  dt  = -c    and   dr  /  dt  =  c  /  L  .  Therefore,  the  propagation  operator  in  (2.32) 
becomes 


d_     ]_d_ 

dz     c  dt      dz 


d      ( _d_     ]_d\ 
dz      L  dr 


V 


/ 


Ldz 


(2.35) 


Hence,  the  full  wave  operator  is 


[dz   cdt 


V2,  +2ik 


\_d_ 
Ldr 


(2.36) 


Introducing  the  dimensionless  transverse  variables  x  =  xv  k  1 2L   and  y  -  yjk  1 2L  ,  so 
that    -^T  =  (k/2L)— -    and    — -  =  (k 1 2L)-^— .   Hence    V\  =  (k / 2L I  Ar  +  — T 

ex2    v      ;55c2         ay2    v      ;c^2  x    v       \ax2    ^2 

d2        d2 
(kl2L)V\ ,  where  V^  = — -  +  -^ .  Substituting  into  (2.36)  and  (2.32)  yields 


dxl     dy' 


4      x     dr 


jr=o. 


(2.37) 


Equation   (2.36)   is   known  as  the   parabolic   wave   equation    [7].   In  dimensionless 
parameters,       the       optical       field       is       represented       by        a  =  \a\e"/>        where 

|a|  =  4k  NeKLE l{y2mc2)  oc  E ,  so  <f  and  a  differ  by  a  multiplicative  constant.  Therefore, 

exchange  <£  for  a,  the  dimensionless  parabolic  wave  equation  (2.37)  becomes 


i 


4      x     dr 


a  =  0. 


(2.38) 


In  the  presence  of  the  electron  beam  current  density,  the  optical  wave  equation 


(2.20)  is 
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V 


2  _  jl_5!_> 

c2  dt2  , 


4x  - 
A(r,t)  = J 


x ' 


(2.39) 


where  J L  is  the  transverse  electron  beam  current  density  [8].  Substituting  (2.30)  into 
(2.39)  yields 


ee 


V\+2ik 


\dz     c  dt  j 


An  - 

J  ±  • 
c 


(2.40) 


dz     c  dt 

Performing    a    dot    product    with    e*  =  (z',1,0)    on   both    sides    of  (2.40)    and    using 
e   ■  e  =  (z,l,0)  •  (-i,l,0)  =  2 ,  yields 


1     2     Jd    1  d^ 

■V;  +  z&  — + 
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ft  k  -      „»    _,a 

j±  •£  e 

c 


(2.41) 


2     A       {dz   cdt. 

The  total  transverse  electron  beam  current  is  the  sum  of  all  single  particle  currents.  The 
single  particle  current  for  one  electron  in  the  beam  is  J±l  =  -ec/3±S{i)(x -Ft)  [8]  where 

fl  is  the  position  of  the  i'h  electron  and  £(3)(x,  -ft)  is  the  three-dimensional  Dirac  delta- 
function.  The  total  beam  current  is 

^i=I^=-«IWP,(x-rl).  (2.42) 

From      (2.10),      the      transverse      motion      is       J3±  =(-  K  /  y\coskoz,sinkoz,0)  = 

We  {(-  K  /  y)iee~,k°:  \.  For  convenience,  we  write   /?±  = iee~ik°2 ,  so  that 

r 


J±-e*e-ia=iecK(e'-e)Y- S{3)(x-r,) 

i         7 


(2.43) 


We  can  simplify  the  sum  TT    by  assuming  a  local  electron  density  p,  so  that  (2.41) 
becomes 


—  V±  +  ik  — + 

2  v  dz    c  dt 


<f  =  -2nieKp\ 


Y 


(2.44) 
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where  C,  =  k0z  +  a     is  the  electron  phase  in  the  combined  undulator  and  optical  fields, 

and  (...)  is  a  weighted  average  of  sampled  electrons  at  a  specific  site  which  replaces  the 

summation  over  all  electrons  within  one  optical  wavelength.  Define  the  dimensionless 
beam  current  j  -  8N(e7r  KL)  p/(y3mc2 ) ,  the  wave  equation  becomes 


4  dz 


a(x>y,z,T)=-j(e-«).  (2.45) 


In  some  cases  where  the  optical  wave  diffraction  is  not  significant,  we  can  ignore  V^  in 
(2.45),  then 

a  =  -j{e-*).  (2.46) 

The  optical  wave  equation  (2.46)  couples  the  electrons  to  the  copropagating  optical  wave 
by  relating  the  electron  phase  to  the  change  in  the  optical  field.  The  equation  (2.46)  can 
be  expressed  in  terms  of  amplitude  and  phase, 

|fl|  =  -y(cos(£  +  0)),  (2.47) 

and       i>  =  /-(sm(^  +  (p)).  (2.48) 

\a\ 

Equation  (2.47)  reveals  that  when  the  electrons  become  bunched  near  (g  +  <f>) «  n ,  the 
optical  field  amplitude  will  grow,  and  when  bunched  near  {C,  +  </>)  «  0 ,  the  optical  field 
amplitude  will  decrease.  Equation  (2.48)  describes  the  optical  phase  evolution. 
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E.         GAIN  DEGRADATION  DUE  TO  ELECTRON  BEAM  QUALITY 

In  the  design  of  an  FEL,  there  is  often  a  trade-off  between  high  beam  current  and 
good  beam  quality  as  determined  by  the  accelerator.  The  dimensionless  current 
j  oc  IN3A]I2  where  /  is  the  actual  beam  current  ,  N  is  the  number  of  undulator  periods, 
and  X  is  the  optical  wavelength.  Usually  the  optical  wavelength  A.  is  fixed  for  a  given 
application  so  that  j  is  maximized  by  adjusting  current  /  and  the  undulator  periods  N. 
Increasing  /  tends  to  degrade  the  beam  quality  from  the  accelerator,  whereas  increasing  N 
tends  to  increase  the  FEL  sensitivity  to  the  beam  quality  because  of  narrowing  the  gain 
spectrum  bandwidth  [9]. 

Emittance  e  describes  one  aspect  of  electron  beam  quality.  A  real  electron  beam 
has  an  emittance  e  describing  how  electrons  enter  the  undulator  with  a  range  of  angles 
6X  with  respect  to  the  undulator  axis,  and  at  a  range  of  transverse  positions  x0  slightly 
off  axis.  The  electron  beam  may  also  have  a  small  energy  spread  Ay  I  y .  The  energy 
spread  and  emittance  of  the  electron  beam  translate  into  a  spread  in  phase  velocities  in 
phase-space.  The  phase  velocity  spread  due  to  energy  spread  is  given  by  [10] 

Av  =4xN  Ayly,  (2.49) 

and  the  phase  velocity  spread  of  an  electron  beam  due  to  emittance  is  given  by 

Av— -^.(rttf+r^).  (2.50) 

Electrons  enter  the  undulator  with  a  spread  in  phase  velocities  due  to  a  spread  in 
energies,  and  are  randomly  distributed  in  both  the  longitudinal  and  transverse  directions. 
The  typical  distribution  of  electron  phase  velocities  is  assumed  to  have  a  Gaussian  form 
about  initial  phase  velocity  v0 , 

/G(<7)=4f— >  (2-51) 


In  a G 

where  oG  =4kN  Ay  I '  y  is  the  standard  deviation  and  q  is  the  phase  velocity  away  from 
v0  distributed  [9].  The  characteristic  function  FG{r)  of  the  distribution  fG(q)  is  defined 
as 
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FG{T)  =  \fG{q)e^dq. 
Evaluating  the  above  integral  equation  for  a  Gaussian  yields 

FG{T)=e-^2'2. 

If  the  beam  quality  is  perfect  so  that  <jg  -  0 ,  then  FG  (r)  =1 .  However,  when  the  beam  is 

not  perfect,  the  characteristic  function  |FG(r)|  decays  over  the  length  of  the  undulator  in 

a  characteristic  time  on  the  order  of  1/  ac  [9]. 

Poor  beam  quality  reduces  the  single-pass  FEL  gain  by  degrading  electron 
bunching.  As  an  example,  consider  an  FEL  design  with  N  =  25  periods  having  a 
dimensionless  current  density  of  j  =10  and  initial  phase  velocity  of  v0  =  3  operating  with 

an  initial  dimensionless  optical  field  a0  =  1.  In  Figure  2.2,  electron  beam  quality  is 
assumed  to  be  nearly  perfect,  i.e.  <jg  =1 .  Some  electrons  gain  energy  and  move  ahead  of 

the  average  flow,   while  other  electrons  lose  energy  to  the  radiation  field  and  move  back 
behind  the  average  flow  causing  visible  spatial  bunching. 


«*  FEL  Pfca..  Space  Evol 

aG=l        ti=0.000 


Figure  2.2.  Phase  space  evolution  for  a  nearly  perfect  electron  beam  quality. 
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In  Figure  2.3,  the  longitudinal  energy  spread  increases  to    <jg  =6,  and  this 

degradation  in  beam  quality  causes  much  less  bunching  when  comparing  to  Figure  2.2. 
The  poor  beam  quality  is  responsible  for  a  decrease  in  final  gain  to  G/=  0.10534  and  a 
smaller  optical  phase  shift  </> . 


3n/2    0 


Figure  2.3.  Phase  space  evolution  for  low  gain  with  electron  beam  degradation. 


F.        TAPERED  UNDULATOR 

The  characteristics  of  the  FEL  interaction  can  be  altered  by  changing  the 
undulator  properties,  i.e.  undulator  polarization,  wavelength,  or  field  strength,  along  its 
length  [5].  As  an  FEL  reaches  saturation,  the  electrons  have  given  up  enough  energy  so 
that  they  are  no  longer  in  resonance  with  the  optical  wave,  which  decreases  the  gain. 
Tapering  changes  the  undulator's  resonance  condition  along  its  length  so  that  the 
electrons  remain  in  resonance  longer,  hence  increasing  the  overall  gain  in  strong  fields 
[7].  From  the  resonance  condition  (2.12),  in  order  to  extract  energy  out  of  the  electron 
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beam,  i.e.  reducing  y,  and  convert  into  radiation  field  and  still  maintain  resonance 
condition,  we  can  either  taper  the  undulator  wavelength  X0 ,  i.e.  reducing  X0,  or  taper  the 

undulator  parameter  K,  i.e.  reducing  the  undulator  field  strength  B.  However,  tapering  K 
is  easier  since  the  undulator  gap  needs  to  be  increased  to  decrease  B.  Tapering  essentially 
has  the  same  mathematical  result  as  accelerating  the  electrons.  This  artificial  acceleration 
is  represented  by  a  dimensionless  parameter  8 .  A  tapered  undulator  represents  an 
advantage  over  the  untapered  undulator  design  when  the  artificial  acceleration  8  exceeds 
the  deceleration  that  can  be  obtained  without  taper  in  strong  optical  fields.  When  the 
undulator  wavelength  is  decreased,  the  artificial  acceleration  8  is  given  by 
8  %  -27rNAA0 1 X0 .       When       the       undulator       field       strength       is       decreased, 

8  «  -4nNK2AB/B(\  +  K2) .  With  tapering  is  included  in  the  undulator,  the  pendulum 
and  wave  equations  become 

g  =  v  =  8  +  \a\cos(£  +  0),  (2.52) 

and  a  =  —j<  e'^  >, 
where  C,  -  [k0(r')dr'  +  kz -cot    [10].  The  wave  equation  remains  the  same,  but  the 

simple  pendulum  equation  includes  an  additional  constant  torque  due  to  the  artificial 
acceleration  8 . 

The  necessary  and  desirable  criteria  for  tapering  is  \a  >  8  >  4b      >  2k  .  The 

condition  on  the  left  is  for  trapping  electrons  in  the  tapered  phase  space;  the  middle 
condition  is  that  the  tapered  acceleration  exceeds  the  natural  untapered  deceleration;  and 
the  condition  on  the  right  is  that  the  tapered  FEL  works  in  the  strong  field  regime  [10]. 

Electrons  initially  near  the  phase  <^0  ~  0  are  accelerated  by  both  the  strong  optical 

field   |<z|    and  the  artificial  "torque".  The  electrons  initially  near   C,0~n    undergo  a 

cancellation  between  the  optical  field  and  the  torque,  leaving  them  trapped  in  closed 
orbits.  The  taper  is  effective  because  electrons  near  the  phase  for  gain  C,  ^n  are  trapped, 
and  continue  to  contribute  to  the  interaction.  Electrons  near  the  phase  for  absorption 
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C,  « 0  are  taken  away  from  resonance  and  eventually  stop  interacting.  The  imbalance 
leads  to  a  net  gain  [10]. 

Compared  to  the  untapered  undulator,  the  tapered  undulator  is  more  efficient  in 
strong  optical  fields,  but  has  smaller  gain  in  weak  optical  fields.  The  efficiency  of  a 
tapered  undulator  is  estimated  at   r\s  =8 /(Sx  N) .  For  example,  if  8  «  280;r   for  an 

undulator  with  23%  taper  in  wavelength  over  N  =  300  periods,  then  the  efficiency  is 
estimated  at  rjs  ^70/(2N )  which  increases  by  a  factor  of  70  over  the  natural  efficiency 

r]  «1/(2tV)  [10].  Tapering  is  a  design  to  obtain  higher  overall  gain  and  efficiency  in 
strong  optical  fields. 

An  example  of  a  tapered  undulator  in  strong  optical  fields  is  shown  in  Figure  2.4 
with  an  untapered  undulator  shown  in  Figure  2.5.  Both  designs  are  based  on  an  undulator 
length  of  L  =  6  m  corresponding  to  N=  300  periods.  The  dimensionless  current  density  is 

j  =  12000  »  k  in  the  high  gain  regime,  and  initial  field  a0  =1 50  Y>k   in  the  strong  optical 

field  regime.  In  Figure  2.4,  the  tapering  rate  8  =  260/r  turns  on  at  zs  =  0.19  along  the 
undulator  resulting  in  a  final  gain  of  Gf  =  556  and  FEL  extraction  efficiency  of 
rj«12.5%.  Initially  the  gain  increases  exponentially  until   zs  =0.19  where  the  taper 

turns  on  and  saturation  is  reached.  With  tapering,  the  gain  continues  to  increase  slowly 
instead  of  decreasing  after  saturation,  as  shown  in  the  upper  right-hand  portion  of  the 
Figure  2.5.  The  phase  velocity  of  the  untrapped  electrons  near  v  «740  is  shifted  by  the 
tapering  acceleration  8  whereas  the  trapped  electrons  remain  near  resonance  due  to 
strong  optical  fields.  The  lower  right-hand  portion  of  the  Figure  2.4  is  the  optical  phase 
shift,  and  this  phase  shift  A0  *  1  An  over  the  length  of  the  undulator  is  a  characteristic  of 
high  gain. 
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***  FQ  Phase  Space  Evolution  *** 
j =12000        aQ=150      vo=0  N-300 

<yQ=l  5=260tt      ti =0.125   ^=0.19 
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Figure  2.4.  Phase  space  evolution  for  high  gain  in  strong  fields  with  a  tapered  undulator. 
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Figure  2.5.  Phase  space  evolution  for  an  untapered  undulator. 
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III.  SIMULATIONS  OF  THE  REGENERATIVE 
MW  FEL  AMPLIFIER 

A.         MW-CLASS  FEL  CALCULATIONS 

In  order  to  destroy  an  incoming  missile,  the  laser  power  required  on  the  missile 
surface  is  about  10  kW/cm  over  a  w0«  5  cm  radius  spot  for  a  two  second  duration.  The 

extinction  coefficient  due  to  aerosols  at  sea  level  is  a  -0.03  km"1,  and  through  e""2 
describes  the  absorption  loss  at  Ul^im  wavelength  over  a  distance  z.  Let  P(  be  the 
Transmitted  Power  that  must  leave  the  ship  to  destroy  a  missile  at  a  range  R,  then  the 
power  at  the  target  is  Pd  =Pt  e~M  .    If  the  range  of  the  missile  is  5  km,  the  laser  power 

which  must  leave  the  ship  is  then  P,  =  [10  kW/cm2]  [it (5  cm)2]  e+003km  '5km  *  910  kW. 
As  an  example,  consider  an  FEL  with  TV  =  300  periods  using  a  beam  of  100  MeV 
electrons  with  a  peak  current  of  400  A,  then  the  peak  electron  beam  power  P  =  400 
Ax  100  MV %  40  GW.  The  micropulses  generated  by  the  FEL  have  a  length  of  0.001  m 
and  are  separated  by  6  m,  then  the  fraction  D  =  0.001/6  denotes  the  duty  cycle  of  the  FEL. 

The  average  power  generated  within  the  FEL  is  given  by  P=PxD«6.7  MW.  The 
efficiency  necessary  to  destroy  an  ASM  is  r)  =  Pt/P«13%.  Therefore,  the  single-pass 
FEL  efficiency  r\  needed  to  supply  the  Laser  Power  to  destroy  the  incoming  missile  is 
13%.  An  approximation  of  the  natural  single-pass  efficiency  of  an  FEL  is  only 
r|  =  1/2N  =  2%  so  that  a  tapered  undulator  with  high  efficiency  is  needed.  The  tapered 

undulator  described  in  the  previous  section  would  be  more  than  adequate.  This  type  of 
illustration  is  the  essential  driving  factor  in  determining  the  FEL  parameters  for  possible 
MW-class  weapon  systems. 
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B. 


LANL  REGENERATIVE  MW  FEL  AMPLIFIER 


The  US  Navy  is  investigating  the  possibility  of  using  a  MW-class  Free  Electron 
Laser  for  ship  self-defense  against  anti-ship  missiles  (ASM).  A  design  workshop  resulted 
in  two  possible  MW  FELs,  the  oscillator  and  the  regenerative  amplifier.  The  MW 
oscillator  design,  although  more  compact  overall,  presents  a  challenge  in  recirculating  an 
intense  electron  beam  because  of  a  phenomena  known  as  coherent  synchrotron  radiation 
(CSR)  feedback  [11].  The  regenerative  amplifier  design  uses  a  longer  undulator  and 
relies  on  higher  extraction  efficiency  to  achieve  high  average  power  so  that  less  current  is 
needed  from  the  accelerator  and  no  recirculation  is  required.  The  regenerative  amplifier 
design  has  been  proposed  by  LANL  and  Boeing  with  the  analysis  being  performed  at 
NPS.  This  section  explores  the  analysis  of  regenerative  amplifier  design  and  optimizes  its 
efficiency. 

A  simple  schematic  of  the  regenerative  amplifier  design  is  shown  in  Figure  3.1. 


Feedback  Loop 


Photo  Injector 


Electron  Dump 


Figure  3.1.  Conceptual  MW  class  Amplifier  FEL. 


The  amplifier  design  does  not  use  energy  recovery  so  there  is  no  requirement  to 
bend  the  high  energy  beam  which  can  cause  deterioration  of  the  electron  beam  quality. 
However,  without  energy  recovery,  higher  extraction  efficiency  is  required.  One- 
dimensional  computer  simulations  are  used  to  describe  a  single  pass  of  the  optical  wave 
and  to  optimize  the  efficiency  using  a  tapered  undulator.  The  tapered  undulator  represents 
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an  advantage  when  the  electron  phase  acceleration  8  exceeds  the  deceleration  that  can  be 
obtained  without  taper  in  strong  optical  fields.  The  single-pass  FEL  efficiency  7  of  a 
tapered  undulator  is  the  fraction  of  the  electron  beam  energy  converted  to  laser  light.  The 
limit  on  77  is  determined  by  the  maximum  taper  rate  8  that  can  maintain  trapped 
electrons  over  the  number  of  periods  N  in  the  undulator  [6].  Gain  G  is  the  fractional 
change  in  power  of  the  optical  field  in  a  single  pass  through  the  undulator.  In  weak  fields 

without  tapering,  the  FEL  amplifier  has  high  gain  described  by  G(r)  *  (l/9)e*J/2)        * 
[12].   The  gain  is  exponential  in  r  along  the  undulator  with  growth  rate  proportional  to 
7'13,    where     j  =  87V \e  n  KL)2  p  F I '  yl mc2     is    the    dimensionless    current    density, 
p=3  x  109 1(A)/  ec7ir2  is  the  electron  beam  density,  and  F  =  "area  of  electron  beam'V'area 

of  light  beam"  is  the  filling  factor. 

Table  3.1    shows  the  proposed  parameters  of  the   MW   regenerative  amplifier 
developed  at  LANL. 


Parameters 

Proposed  Value 

Beam  energy     E 

100  Me V 

Beam  radius     rh 

0.17  mm 

Energy  spread    Ay  1  y 

0.02% 

Pulse  duration 

20  ps 

Pulse  frequency 

500  MHz 

Peak  current     / 

400  A 

Average  current  / 

0.2  A 

Undulator  parameter  K 

1.71 

Undulator  length  L 

6  m 

Number  of  undulator  periods  ./V 

300 

Undulator  period  A,n 

2.0  cm 

Optical  wavelength  X 

1.0  urn 

Table  3.1 .  Parameters  for  the  MW  regenerative  amplifier  design  developed  at  LANL 
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1.  ID  Efficiency  Optimization 

The  proposed  M W  RAFEL  uses  a  beam  of  1 00  MeV  electrons  with  a  peak  current 
of  400  A,  yielding  a  peak  electron  beam  power  of  40  GW.  The  average  current  of  the 
RAFEL  is  0.2  A,  so  that  the  average  electron  beam  power  is  20  MW.  The  single-pass 
FEL  efficiency  t]  is  the  fraction  of  power  extracted  from  electron  beam  in  one  pass 
through  the  undulator.  The  RAFEL  will  feedback  on  the  order  of  0.01%  to  1%  of  the 
optical  power  and  require  an  extraction  efficiency  of  approximately  10%  to  15%  to 
provide  2  MW  to  3  MW  of  optical  power  in  the  infrared  (IR). 

In  an  attempt  to  achieve  the  desired  efficiency,  the  initial  undulator  length  is  L  = 
4m  corresponding  to  N  =  200  periods.  The  dimensionless  initial  optical  field  a0  and 

tapered  undulator  of  strength  S  are  varied  to  find  the  optimum  efficiency.  The  maximum 
extraction  efficiency  found  was  8%  which  is  less  than  the  desired  efficiency. 

In  an  effort  to  improve  the  design  efficiency,  the  length  of  the  undulator  is 
increased  from  4  m  to  6  m.  The  change  in  undulator  length  corresponds  to  dimensionless 
current  density  of  j  «  12,000  with  N  =  300  undulator  periods.  The  values  of  a0  and  S 
are  varied  to  find  the  optimum  efficiency  using  numerical  simulations.  The  electron 
beam  has  a  Gaussian  spread  crG  =1.0,  and  the  field  strength  is  a0  =  1 50  with  taper  rate 

5  =260;r  starting  at  time  ts  =0.19  along  the  undulator.  Figure  3.2  is  the  3D  plot  of 
efficiency  versus  a0  and  8  .  The  maximum  efficiency  rjmax  is  13.3%  with  gain  G  =  556. 
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Efficien  cy  Optimization 


300 


0     0 


Figure  3.2.  3D  plot  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Tapered  Field  Strength. 
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The  effect  of  changing  a0  away  from  its  optimum  value  is  smaller  than  changing  8  away 
from  its  optimum  value  as  shown  in  Figure  3.3. 


400 


Efficiency  Optimization 


300 


Fig  3.3.  Contour  plot  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Tapered  Field  Strength. 
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For  a  better  understanding  of  the  physical  mechanism  underlying  the  operation  of 
this  MW  RAFEL,  examine  Figure  3.3  closely  in  five  regions.  Region  I  is  the  "optimum 
region".  As  shown  in  Figure  3.4,  the  taper  turns  on  at  zs  =0.19,  and  approximately  half 

of  electrons  were  trapped  near  resonance  v0  =  0 ,  resulting  in  an  acceptable  efficiency  of 

77  =13.3% .       The       final       gain       Gf  =  556       corresponds       to       feedback       of 

/  =  1  l(Gj  + 1)  =  1 .8  x  1 0"3 .  The  large  optical  phase  shift  <f>  in  the  lower  right-hand  portion 

of  Figure  3.4  represents  good  optical  guiding  which  is  an  important  part  of  a  high  gain, 
high  current  FEL. 
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Figure  3.4.  Phase  space  evolution  for  the  MW  RAFEL  at  optimum  efficiency. 
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Region  II  is  the  "high  tapering  rate  region".  Increasing  the  tapering  rate  8  beyond 
the  optimum  value,  the  efficiency  drops  to  8.3%  due  to  fewer  electrons  being  trapped  near 
resonance  v0  =  0  as  shown  in  Figure  3.5.  The  taper  turns  on  at  rs  =0.19  as  before.  The 
final  gain  Gf  =  347  corresponds  to  feedback  of  /  =  1/(G/  +I)  =  2.9xl0  3.  <fi  is 
relatively  small  which  would  cause  less  optical  guiding. 
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Figure  3.5.  Phase  space  evolution  for  the  MW  RAFEL  with  high  tapering  rate. 
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Region  III  is  the  "low  tapering  rate  region".  For  S  below  the  optimum  value,  the 
efficiency  drops  to  6.5%  due  to  insufficient  taper  rate  as  shown  in  Figure  3.6.  While  many 
electrons  are  trapped  near  resonance  v0  =  0 ,  the  low  taper  rate  does  not  extract 
significant  energy  from  them.  The  taper  again  turns  on  at  zs  =0.19.  The  final  gain 
G,  =294  corresponding  to  feedback  of  f  =  \/(G,  +1)  =  3.4x1  (T3 . 
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Figure  3.6.  Phase  space  evolution  for  the  MW  RAFEL  with  low  tapering  rate. 
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Region  IV  is  the  "high  optical  field  region".  Starting  with  high  optical  field 
strength  a0 ,  the  separatrix  does  not  have  a  chance  to  grow  and  trap  electrons,  resulting  in 
a  lower  efficiency  of  9.5%  as  shown  in  Figure  3.7.  The  taper  turns  on  at  rs  -  0.16 .  The 
final  gain  Gf  =144  corresponds  to  feedback  of  f  =  \l{Gf  +  l)  =  6.9x  10"3 . 
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Figure  3.7.  Phase  space  evolution  for  the  MW  RAFEL  with  high  optical  field  strength. 
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Region  V  is  the  "low  optical  field  region".  Starting  with  low  optical  field  strength 
a0 ,  the  initial  separatrix  is  too  small  so  that  more  electrons  leak  out  from  the  separatrix, 
resulting  in  a  lower  efficiency  of  10.9%  as  shown  in  Figure  3.8.  The  taper  turns  on  at 
r5=0.32.  The  final  gain  Gf  =25,700  corresponds  to  feedback  of 
/  =  1/(G/+I)  =  3.9xl0"5. 
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Figure  3.8.  Phase  space  evolution  for  the  MW  RAFEL  with  low  optical  field  strength. 

The  above  results  for  efficiency  should  be  considered  an  over-estimate  since  the  effect  of 
diffraction  of  the  optical  beam  or  longitudinal  pulse  effects  are  not  included. 
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2.  3D  Diffraction  Optimization 

From  Chapter  II,  the  parabolic  wave  equation  can  be  written  as 


4      i     dr 


a(x,y,z,T)=-j(e-*}.  (3.1) 


Equation  (3.1)  governs  the  dynamics  of  the  optical  wave  over  many  optical  wavelengths  in 
the  longitudinal  z-dimension,  and  the  V^  operator  properly  describes  the  diffraction  of  the 
optical  wave  in  the  transverse  (x,y)  directions.  A  coherent,  freely-propagating  optical 
wave  of  radius  w0  will  spread  due  to  natural  diffraction  as  it  traverses  the  undulator  length 
L.  The  Rayleigh  length,   nw2oIX,  is  the  characteristic  distance  over  which  the  initial 

optical  mode  area  doubles  in  size.  If  the  optical  mode  area  becomes  much  larger  than  the 
electron  beam,  the  optical  field  will  decouple  from  the  electron  beam.  Therefore,  it  is 
important  that  the  Rayleigh  length  is  comparable  to  the  undulator  length  L,  and  the  optical 
mode  waist  w0  is  not  too  much  larger  than  the  electron  beam  radius  r&  [13]. 

Three-dimensional  simulations  are  used  to  describe  a  single  pass  of  the  optical 
wave  through  the  undulator.  Initially,  the  optimal  values  from  one-dimensional  simulation 
were  used  as  the  inputs  for  a  3D  simulation.  The  resulting  3D  simulation  is  shown  in 
Figure  3 .9  with   ao  =  1 50 ,  8  =  260;r ,  j '  =  2.6  x  104 ,  and  the  taper  turning  on  at  rs  =  0.2 . 

The  evolution  of  the  optical  mode,  |a(x,r)|  shown  in  the  upper-left  window  of  Figure  3.9, 

is  dominated  by  diffraction  after  the  taper  is  turned  on  at  rs  =  0.2 .  The  intense  electron 

beam  provides  guiding  over  the  first  one  fifth  of  the  undulator,  but  after  tapering  begins 
diffraction  overcomes  this  focusing  over  the  rest  of  the  undulator.  The  top-center  window 
shows  a  cross-sectional  view  of  the  optical  mode,  |a(x,_y)|,  at  the  end  of  the  undulator, 

and  this  optical  mode  has  a  drop  in  optical  power  in  the  center.  This  drop  is  due  to  a  loss 
of  guiding  of  the  optical  mode  [14].  The  window  in  the  top-right  lists  the  dimensionless 
parameters  for  this  simulation;  ax  and  a    correspond  to  the  radial  size  of  the  electron 

beam,  while  <Je  and  aG  account  for  beam  quality. 
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Figure  3.9.  3D  diffraction  simulation. 

The  graph  at  the  right-middle  shows  the  development  of  the  optical  phase  along  the 
undulator  length.   At    r,  =0.2,  the  optical  phase  change  along  the  undulator  length 

becomes  negative  and  drops  rapidly  which  causes  the  optical  wavefront  to  be  excluded 
from  the  electron  beam  decreasing  the  interaction  strength.  The  lower-right  window  is  the 
plot  of  the  dimensionless  optical  power  growth  along  the  undulator,  and  the  plot  of  the 
natural  logarithm  of  the  single-pass  gain  along  the  undulator.  At  the  middle-left  shows  the 
bunching  current  in  the  electron  beam  development  along  the  undulator,  with  the  end  view 
in  the  center.  The  bottom-left  plot  is  the  electron  beam  phase  velocity  evolution  along  the 
undulator,  and  at  the  bottom  center  shows  the  electron  phase  space  plot  at  the  end  of  the 
undulator.  Initially  the  optical  guiding  confines  the  optical  mode  until  the  taper  turns  on  at 
rs  =  0.2,  then  diffraction  effects  begin  to  spread  the  beam.  The  optical  phase  begins  to 

grow  until  the  taper  turns  on,  then  the  change  in  optical  phase  is  negative  as  the  guiding 
effect  is  lost.  The  bottom-center  window  shows  that  at  the  end  of  undulator,  no  electrons 
are  trapped  at  resonance,  and  the  tapering  effect  has  failed. 
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In  order  to  search  for  an  optimum  efficiency  using  numerical  3D  simulations,  the 
initial  field  a0,  the  taper  rate  8 ,  and  the  time  to  turn  on  the  taper  ts  are  varied.  The 
optical  field  amplitude  a0  is  varied  from  50  to  275  in  steps  of  25.  For  each  value  of  a0,  the 
taper  field  rate  8  is  varied  from  30;r  to  330^  in  steps  of  25^-,  and  for  each  value  of  a0 
and  8 ,  the  turn-on  time  rs  is  varied  from  0. 1  to  0.5  in  steps  of  0.02.  First  we  explore  the 
MW  RAFEL  efficiency  by  varying  8  and  r.  at  fixed  values  of  a0. 


3D  Efficiency  Optimization  (for  a0  ~  50) 
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Figure  3.10.  3D  and  Contour  plots  of  Efficiency  vs.  Tapered  Field  Strength 
and  Taper  turn-on  time  at  a0  =  50. 
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Figures  3.10  through  3.14  show  the  dependence  of  efficiency  rj  as  rs  and 

8  are  varied.  These  two  parameters  can  have  a  significant  impact  on  efficiency 
when  moved  away  from  their  optimum  values.  The  two  graphs  of  each  figure  are 
plotted  with  the  same  data  with  the  top  one  is  in  3D  and  the  bottom  one  is  in  2D. 


3D  Efficiency  Optimization  (for  ao=100) 
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Figure  3.1 1.  3D  and  Contour  plots  of  Efficiency  vs.  Tapered  Field  Strength 
and  Taper  turn-on  time  at  a0  =  100. 
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At  the  initial  optical  field  amplitude  of  a0  =  50,  the  efficiency  peaks  at  rjmm  m  6.1% 
at  values  of  8  =  200/r  and  vs  =  0.28 .  As  a0  increases,  the  efficiency  peaks  at  slightly 
higher  efficiency  rjmax  «7.0%  at  lower  values  of  5  =  \90n:  and  ts  =0.21  as  shown  in 
Figure  3.12. 


10-r 


tK%)5 


3D  Efficiency  Optimization  (for  a0=175) 
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Figure  3.12.  3D  and  Contour  plots  of  Efficiency  vs.  Tapered  Field  Strength 
and  Taper  turn-on  time  at  a0  -  175. 
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The  peak  efficiencies  ?jmax  continue  to  increase  to  a  optimum  efficiency 
rime*  w  1-2%  corresponding  to  optimum  values  a0  =  225,  8  =  \S0tt  ,  and  taper  turn-on 
time  of  ts  =0.18  as  shown  in  Figures  3.13. 


3D  Efficiency  Optimization  (for  »  =225) 
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Figure  3.13.  3D  and  Contour  plots  of  Efficiency  vs.  Tapered  Field  Strength 
and  Taper  turn-on  time  at  a0  =  225. 
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Increasing  optical  field  strength  a0  beyond  its  optimum  value,  the  separatrix  does 
not  have  a  chance  to  grow  and  trap  electrons,  resulting  in  a  slightly  lower  efficiency  of 
7.1%  as  shown  in  Figure  3.14. 


3D  Efficiency  Optimization  (for  ao=250) 
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Figure  3.14.  3D  and  Contour  plots  of  Efficiency  vs.  Tapered  Field  Strength 
and  Taper  turn-on  time  at  a0  =  250. 
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Next,  explore  the  effect  on  efficiency  by  varying  a0  and  zs  at  a  fixed  value  of  8  . 
As  shown  in  Figure  3.15,  starting  at  a  low  value  of  the  taper  field  rate  £  =  55/r,  the 
maximum  efficiency  obtained  was  3.2%  with  ri  »0.2.  This  low  efficiency  is  due  to 
insufficient  taper  rate. 


3D  Efficiency  Optimization  (for  8=557C) 
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Figure  3.15.  3D  and  Contour  plots  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Taper  turn-on  time  at  5  =  55tt  . 
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As  8  increases,  the  efficiency  increases  to  t]max  »  5.3%  at  a  slightly  lower  ts  as 
shown  in  Figure  3.16.  The  peak  efficiency  increases  to  7max  «7.2%  corresponding  to 
optimum  values  of  8  =  1 807T ,  a0  =225,  and  an  earlier  taper  turn-on  time  of  r  =0.18. 
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3D  Efficiency  Optimization  (for  5=^1057t) 
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Figure  3.16.  3D  and  Contour  plots  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Taper  turn-on  time  at  8  =  105/r . 
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As  shown  in  Figure  3.17,  when  the  optimum  value  of  taper  rate  is  used,  8  =  180tt  , 
varying  the  field  a0  has  little  effect  on  efficiency,  whereas  changing  the  taper  turn-on  time 
from  its  optimum  value,  xs  «18,  results  in  a  significant  efficiency  drop. 


3D  Efficiency  Optimization  (for  6=180rc) 
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Figure  3.17.  3D  and  Contour  plots  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Taper  turn-on  time  at  8  =  1 80;r . 
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Increasing  8  beyond  its  optimum  value  causes  the  efficiency  to  drop  significantly 
due  to  fewer  electron  being  trapped.  As  shown  in  Figure  3.18,  the  efficiency  drops  to 
only  1.9%  at  r  =  0.22  with  8  =  280;r . 


3D  Efficiency  Optimization  (for  6~2807t) 
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Figure  3.18.  3D  and  Contour  plots  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Taper  turn-on  time  at  8  -  280;t  . 
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Now  consider  the  three-dimensional  simulation  of  a  single  pass  of  the  optical  wave 
through  the  undulator  at  the  optimum  values  of  8  =  180;r,  zs  =0.18,  and  a0  =  225.  As 

shown  in  Figure  3.19,  optical  guiding  occurs  and  stays  with  the  intense  electron  beam 
through  the  undulator  as  seen  in  upper-left  and  middle-left  windows.  The  top-center 
window  shows  a  much  tighter  cross-sectional  end  view  of  the  focused  optical  power  than 
seen  in  Figure  3.9.  The  bottom-center  window  shows  approximately  50%  of  the  electrons 
are  trapped,  resulting  in  optimized  gain  and  efficiency.  The  middle-right  window  shows 
the  optical  phase  initially  increasing  until  the  taper  turns  on  at  zs  «  0.18 ,  then  decreasing 
slowly  The  FEL  performance  in  Figure  3.19  is  much  improved  over  that  of  Figure  3.9.  In 
steady-state    operation,     the    final    gain     G/=616     corresponds    to    feedback    of 

/  =  1/(G/+I)  =  1.6xl0-3. 
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Figure  3.19.  3D  diffraction  simulation  at  optimum  values 
5  =  \S0tt,  t.  =  0.18,  and  a0  =  225. 
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IV.      CONCLUSIONS 

The  development  of  a  speed-of-light  hard-kill  weapon  system  for  military 
applications  may  represent  a  significant  advancement  in  technology  over  present 
conventional  kinetic  weapon  systems.  The  FEL  has  the  potential  for  use  as  a  weapon 
system  against  anti-ship  and  theatre  ballistic  missiles.  The  infinite  magazine,  rapid 
response,  and  wavelength  tunability  at  high  powers  make  it  suitable  and  desirable  for 
shipboard  self-defense  and  for  land-based  theater  ballistic  missile  defense. 

Three  dimensional  computer  simulations  indicate  that  the  proposed  LANL 
Regenerative  MW  FEL  design  does  not  yield  the  desired  extracted  efficiency  of  rj  «  1 5%  . 
The  optical  field  amplitude  a0,  taper  rate  5 ,  and  taper  turn-on  time  rs  were  varied  to 
search  for  the  optimum  efficiency,  resulting  in  maximum  efficiency  of  only  rjmax  «  7.2% . 
Changing  S  or  zs  away  from  their  optimum  values  has  a  significant  effects  on  the  MW 

RAFEL  design  efficiency,  but  changes  in  a0  do  not  have  much  of  an  effect.  The  efficiencies 
found  here  are  still  an  over-estimate  since  they  do  not  include  longitudinal  pulse  effects.  It 
remains  to  be  determined  whether  a  smaller  efficiency  is  acceptable  or  if  1 5%  efficiency  is 
actually  necessary  in  the  FEL  weapon  design.  If  15%  efficiency  is  necessary,  a  longer 
undulator,  perhaps  8  or  10  meters  may  be  sufficient  to  reach  that  goal. 
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